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Rich experimental data accumulated in the past few years in the hadronic decays allow one 
to check the quark combinatorics relations for a new type of processes, namely: quark jets in the 
decays Z" qq ^ hadrons. In this paper we review quark combinatorics rules for the yields of 
vector and pseudoscalar mesons, V/ P, in the central and fragmentation regions of the hadronic 
decays. It is stressed that in the central region a direct verification of quark combinatorics rules 
, is rather problematic because of a considerable background related to the decay of highly excited 

' resonances; however such a verification is possible in the fragmentation region, at Xhadron ^ 0.5 — 1, 

I where the contribution of resonance decays is suppressed due to the rapid decrease of spectra with 

increasing Xhadron- It is shown that the experimental data in the fragmentation region for p^/tt" and 
p/n'^ are in a reasonable agreement with the predictions of quark combinatorics. The ratios of the 
heavy meson yields, B* / B and D* /D, are also discussed: the data demonstrate a good agreement 
with the quark combinatorial results. We analyse the structure of the suppression parameters for 
strange and heavy quark productions in soft processes and estimate their order of magnitudes in 
^SJ , the multiperipheral processes; the ratio i^^/yr* at xhadron ~ 0.5 — 0.8 and production probabilities 

of the cc and bb mesons are in a good agreement with the estimates. 
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I. INTRODUCTION 



00 I The relations of quark combinatorics for hadron yields were originally based on the qualitative treatment of a 
multiparticle production process as a process of production of a "cloud" of non-correlated quarks and antiquarks 
followed by their occasional fusion into mesons and baryons In this picture the quark combinatorics approach 

, echoes the statistical treatment of reactions of nucleus decays Q, and actually it was initiated by the success of 
' the statistical methods in the description of multinucleon processes. There is one more common point in treating 
Oh multinucleon and multiquark processes: the stage of final state interactions. In multinucleon reactions the Watson- 
Migdal factor Q leads to an increase of the low-energy contributions in the nucleon-nucleon final state interactions. 
In multiquark reactions the analogous effect of the qq and qqq final state interactions, as one can guess, should enlarge 
the production of the comparatively low-mass mesons and baryons: in |^] a hypothesis was made about the dominant 
^ ' contribution of several low-mass qq and qqq multiplets to the spectra of produced hadrons. 

This particular story reflects a general trend of how methods and approaches used in nuclear physics influenced the 
physics of strong interactions. But, as time goes on, our understanding of QCD, the theory of strong interactions, 
becomes deeper, and now a new level of understanding of quark-gluon physics in the soft interaction region forces us 
to review the relations of quark combinatorics. 

a) Quark combinatorics for hadron production. 

The rules of quark combinatorics, as it was stressed above, were initially suggested for finding probabilities of 
hadron yields in the multiparticle production processes Later, they were extended to meson decay processes, 

see Q . In ||,^ (see also ) the hadronic production processes were treated on the basis of the following qualitative 
picture. In the multiparticle production process there is a stage of a cloud of constituent quarks. When joining each 
other, they create mesons (qq), baryons (qqq) and antibaryons (qqq). The probabilities to form hadrons of different 
sorts are defined by combinatorial rules which have been obtained under the assumption that quarks and antiquarks 
of the sea are produced uncorrelated and there is neither spin nor isospin alignment. 

Our knowledge of the multiparticle production mechanism, i.e. the structure of amplitudes responsible for the 
inclusive production, is now, however, more complete, thus allowing us to investigate the basis of quark combinatorics 
at a new level. In the present paper, we investigate 

(i) the ratio of vector/pseudoscalar meson yields both for the central and fragmentation regions of the quark jets in 
the hadronic Z° decays (Sections 2 and 3, respectively), and 

(ii) the structure of the suppression parameters for production of strange and heavy quarks in the multiperipheral 
processes (Section 4). 
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Time is ripe for this review, for in the meantime quark combinatorics has become widely used for the consideration 
of meson resonances with masses 1.0-2.5 GeV, with a goal to determine the quark-gluon content of these states by 
investigating the transitions resonance mesons In the investigation of exclusive decay processes [0-^, 

the qualitative picture of the cloud structure of sea quarks is not used (that makes the conclusions more model- 
independent). However, the notion "suppression parameter" for the production of strange quarks plays an important 
role , and it is the same as what is used in multiparticlc production processes. Thus, we have to re-think quark 
combinatorics on a new basis both theoretically and experimentally. 

b) Hadronic Z° decays. 

We discuss here the yields of vector {p^,uj, K*) and pseudoscalar (tt, K) mesons in the hadronic decay. Currently 
there exists rich experimental information on these processes, which are determined by transitions ^ qq ^ 
hadrons. First, we consider the production of light-flavour mesons, which are created in the quark jets Z° uu, 
Z° dd, Z" ss with probabilities 0: 

uu:dd:ss~ 0.26 : 037 : 0.37 (1) 

The large mass of the boson makes it possible to observe in hadronic decays Z'^ ~+ qq —+ hadrons the characteristic 
features of both multiparticlc production (central region of quark jets) and meson decay processes (fragmentation 
production). 

The data accumulated by the Collaborations ALEPH L3 H, DELPHI (Ts) and OPAL Q provide spectra of 
vector and pseudoscalar mesons, da/dx{Z^ — > F -f X) and da/dx{Z^ ^ P + X) in a broad interval of x, thus giving 
an opportunity to compare them with the predictions of quark combinatorics. 

c) Hadron production in the central region, and the vector-pseudoscalar ratio. 

In Section 2 the hadron production in the central region in the quark jet is considered. We discuss the prompt 
vector meson (V) and pseudoscalar meson (P) yields . For particles belonging to the same qq multiplet we obtain: 

Pprompt 

This is a well-known result of the quark combinatorics 0]^,^ for hadron-hadron collisions. Hence, our analysis shows 
that the ratio Vprompt / Pprompt is the same for pomeron ladder in hadron-hadron collisions and for quark jets, despite 
of the different structure of the colour exchanges in these processes. 

As we have already said, the rules of quark combinatorics for hadron-hadron collisions were introduced in ||l|,§,D. 
Investigations of the QCD-pomeron shed light on quark-gluon structure of the multiperipheral ladder in hadron 

collisions and allowed us to deal with meson yields in the central region on a new level. Therefore, in Appendix A we 
give a new analysis of the results of quark combinatorics for Vprompt / Pprompt at x '--^ 0, using Lipatov's pomeron 
as a guide. 

d) Suppression of strange quark production. 

An important feature in hadron production is the suppression of the production probability of the strange quark 
in the multiperipheral ladder: 

uu : dd : ss = 1 : 1 : X , (3) 

with < A < 1. This topic is discussed at Section 4: it is shown that A ~ m'^/m'^, where m and are masses of 
non-strange and strange constituent quarks, respectively. 

This suppression results in a suppression of the production of mesons containing strange quarks that is observed in 
the hadronic Z° decays. 

e) Soft colour neutralisation of quarks in the hadronic decay. 

When considering central production at the decay Z° — s- — s- hadrons we start with the standard mechanism of 
soft colour neutralisation of the outgoing quarks: newly-born quark-antiquark pairs are produced in multiperipheral 
ladder (see Fig. la), which provides the transfer of the colour from antiquark to quark. The discontinuity of the 
self-energy diagram of Fig. lb (determined by cutting through hadronic states, dashed line in Fig. lb) defines 
the transition cross section Z'^ —^ hadrons, while the quark-gluon block inside the big quark loop determines the 
confinement forces. 

Likewise, the inclusive production cross section of the meson in the central region is provided by the discontinuity 
of the diagram of Fig. Ic. The quark loop meson qq ^ meson shown in Fig. Ic with the production of vector 
or pseudovector mesons determines relative probabilities of these particles. The chain of the quark loops shown on 
Figs, lb, Ic and Id (below we denote this chain as A) contains both colour singlet (c = 1) and colour octet (c — 8) 
components: A = Ai + Ag. According to the rules of 1/A^-expansion the main contribution is due to the octet 
component. The idea that the quark leaves the confinement trap by a production of new quark-antiquark pairs is 
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rather old; it has been discussed long ago (for example, see ||^, Sections 7 and 9, as well as [Q). A new impulse in 
understanding the confinement mechanism has been brought by Gribov jl^]. Following the ideas of we use the 
jet structure shown in Fig. la assuming that the t-channel exchange by quark is a constructive element of the jet. 

Spectroscopic calculations (see, for example, po| ) support the hypothesis about the scalar type of confinement 
forces; accordingly, we assume that the chain A realises the i-channel exchange with ~ 0+. 

The calculation of the block for central meson production performed in Section 2 proves that equation (2) is fulfilled, 
if the wave functions of mesons {V and P) belonging to the same multiplet are equal. 

f ) Prompt hadron production and production due to the decay of high resonances. 

The equality (2) is valid for prompt meson production, while the decays of highly excited states violate this ratio 
as is seen in the experiment, providing the increase of the rate of light mesons due to resonance decay. As to p/n and 
K*/K, the decays increase the contribution of pseudoscalar component. According to /n^ — 0.15 ± 0.03, and 

K* /K — 0.40 ± 0.06, that indicates on a large contribution into spectra from the decays of highly excited states. 

The problem of the production and decay of highly excited states in hadron-hadron collisions had been discussed 
in pl[ |. The conclusion was similar: average multiplicities of the produced light mesons and baryons result mainly 
from the cascade decays of the highly excited resonances. 

The multiperipheral structure of the ladder allows us to estimate the resonance masses which are important for 
meson production; they are of the order or less than 2 GeV (Section 5) 

The existence of the decay of highly excited resonances is a reality that one should take into account in the 
verification of quark combinatorial rules. We discuss several ways of solving this problem. 

g) Heavy meson production in — + 66 and Z° — + cc decays. 

One way is to check quark combinatorics for heavy particle yields, where the cascade multiplication is suppressed. 
An ideal example could be the production of mesons containing a 6-quark. In fact, for the beauty mesons the ratio 
B* / B observed in the experiment agrees with (2). When the lowest 5- wave multiplet dominates in the production 
of heavy mesons, then one has (provided the equation (2) is fulfilled): B ~ Bprompt + B*^^^^^ = ABprompt, and the 
ratio of vector and pseudoscalar mesons is ~ 0.75. 

At the experiment on Z° decay it was observed: B* / B ^ 0.771 ± 0.075 H, 0.72 ± 0.06 H, 0.76 ± 0.10 [H, 
0.76 ± 0.09 [ll; the mean value is 0.75 ± 0.04. 

For charmed mesons D* / D = 0.60 ±0.05 0.62 ±0.03 |3, 0.57 ±0.05 |||. The mean value is 0.61 ± 0.03,which 
means a rise of the contribution from the decay of the non-5-wave multiplets. 

One should stress that the ratios V/ P for beauty and charmed mesons are saturated in the fragmentation region 
due to the transitions — s- 66 and — > cc. Therefore, in Section 3 we re-analyse quark combinatorics for the 
fragmentation region of the hadronic Z^ decay. 

The production cross section of mesons in the fragmentation region is determined by the discontinuity of the 
diagram of Fig. Ic (the cutting of ladder diagram is shown by dashed line). Direct calculations demonstate that (2) 
is fulfilled with rather good accuracy for the fragmentation region as well, provided the wave functions of vector and 
pseudoscalar mesons are equal. 

h) Production of Hght-flavour mesons in the fragmentation region of the hadronic decay. 
Investigations of meson production in the fragmentation region open the way to test the rules of quark combinatorics 

for light-flavour hadrons, and to verify (2) in particular. As is said above, in the spectra of light -flavour hadrons the 
contribution of the component related to the decay of highly excited states dominates. Still, in the case of jet processes 
this component dominates in the central region, at a; ~ 0, but not in the fragmentation one. The hadronic spectra for 
jets are maximal at a; ~ 0, and they decrease rapidly with the growth of x. As a result, the component which comes 
from a resonance decay decreases quickly, because the decay products share the value Xresonance, thus entering the 
region of smaller x. In due course this effects a fast growth of relative contributions from prompt particle production. 
Therefore, the measurements of particle yields at a; '--^ 0.5 — 1.0 provide an opportunity for model-independent testing 
of quark combinatorics. 

In Section 3 we compare the ratios p^/vr" and K* / K measured in [Tl| ] at large x with quark combinatorics predic- 
tions; in Section 5 the baryon-meson ratio, -pjix^ ^ is discussed, comparison demonstates reasonable agreement of the 
experimental data with the predictions of quark combinatorics. 

i) Heavy quark production suppression. 

Multiperipheral dynamics of the quark-gluon ladder in the processes of Fig. 1 results in a strong suppression 
of production of new heavy quark pairs, QQ. The suppression parameter (relative probability) is of the order of 
Aq ~ vr? I (toq ln^(A^/r7iQ)), where mq is the mass of heavy quark and A is the QCD scale parameter, A ^ 200 MeV. 
One has Ac ~ 2.8 • 10"^ ^nd Ab ~ 1.1 • lO'^. 

The inclusive production of the cc and 66 mesons in decays agrees with this estimation (Section 4) . 

j) The ratio Loj p^ and interference of flavour amplitudes at a; ~ 1. 



3 



The rules of quark combinatorics give us ^^prompt / Pprompt = 1- It can be proven that the decays of highly excited 
states effect weakly this equality. Indeed, the masses of p and uj are almost equal, so the phase spaces of the decay 
processes influence equally the probabilities of their yields. The experiment confirms this statement: at x ^ 0.1 — 0.2 
the value lu/p'^ = 1 is observed. However, at larger x the experimental value is uj/p^ < 1. This reveals the interference 
of flavour amplitudes. The fact is that the decay vertices for Z'^ uu and Z° dd have opposite signs, thus 
enhancing the coherent production of the p° meson and suppressing the production of uj (Section 6). 



II. PROMPT PRODUCTION IN THE CENTRAL REGION OF THE QUARK JET: THE RATIO V/P = 3 

In this Section we discuss the prompt production of vector and pseudoscalar mesons (for definiteness, p and tt) 
in the central region of a quark jet. The production cross section is determined by the discontinuity of the diagram 
shown in Fig. Ic; it is re-drawn in Fig. 2a. A particular feature of the production of p and tt is the presence of a loop 
diagram, which is shown separately by Fig. 2b. Below we calculate these loop diagrams for p and tt using the spectral 
integration technique which is discussed in detail in p9| , ^ . Within this technique the loop diagrams are expressed 
in terms of the p and tt light-cone wave functions. 

But first, let us present the result of our calculations. 

Direct calculations lead to the following formulae for the inclusive cross section of the p and tt mesons at x ^ 0: 



1 

-> - + ^) = / / rf^k,v4(e, k,) . u.iw?, wi) . (4) 



Here ijjp and V'tt are quark wave functions of p and tt mesons, ^ and kj^ are quark light-cone variables (the fraction 
of momentum carried by a quark along the z-axis and its momentum in the (x,y)-plane, respectively). In (4) one 
can see explicit expressions related to the production of p and tt mesons. The rest (contribution from the large quark 
loop as well as from ladder diagrams) is designated in (4) as Ilz{Wi ,W2), which depends on the invariant energies 
squared for the quark chains, Wi and W|. Multiperiphcral kinematics give: 

W^W^^a^-0{m^ + k'^)Ml. (5) 

Here Mz is the mass of the boson. 

The factor 3 in the p production cross section is the result of summing over polarizations of the vector particle. 

Equation (4) directly demonstrates that, if quark wave functions of p and tt arc identical, that is assumed by the 
quark multiplet classification of these mesons, then at x ~ wc have da{Z^ ^ P + X)/dx : d(j{Z^ tt + X)/dx — 
3:1. Let us stress again: the equation (4) as well as diagrams of Figs. Ic, 2a, 2b stand for promptly produced 
mesons. 



A. Spectral representation for the loop diagram of Fig. 2a 

The calculation of the diagram of Fig. 2a in terms of spectral integration is performed in the following steps (see 
also glo)): 

(i) quark loops in Fig. 2a are taken energy-off-shell; 

(ii) for the energy-off-shcU quark loops the discontinuities are calculated (corresponding cuts arc shown by the dotted 
lines I, II, HI and IV); 

(iii) the spectral integrals are determined by the discontinuities being the integrands. 

First, consider the double spectral integral which corresponds to the cuts III and IV, which are the spectral integrals 
over effective masses squared, NP and Af ^, in the transitions Z^ — > qq and qq ^ Z^: 

oo 



dMMM'^ gz{M^)9z{M'') 



j 7r2 (M2-M|)(Af'2-Af|) 

4m^ 



2— —2 -2^ d^{Pz;Pl,P2)d^{Pz;P3,P4) 



X Sz Tip; q^,q2)A{Wt,qi)A{Wi,qi) . (6) 
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Here gziM"^) is the vertex function for the qq transition and Mz is the Z'^ boson mass. The factors d^{Pz;pi,P2) 

and d^{P'^;p3,p4) are phase spaces related to the cuts III and IV. Sz is the spin factor for the big quark loop in Fig. 
2a. The amplitudes A{Wf , qf), {i = 1, 2), refer to the quark-gluon chains with the i-channel scalar quantum numbers 
(see discussion in Section l.e), and T{p;qi,q2) is the block which corresponds to the small quark loop (the qq loop 
for production of p and tt mesons). 

The characteristic feature of the spectral integral (6) is the large value of Mz- Because of that one can replace, 
with rather good accuracy, the poles of the spectral integrand by half-residues: 

z^<5(A/2-M|), —J—-^^nSiM''~Ml). (7) 



Equation (6) stands for the discontinuity of the amplitude that results in different signs for the half- residues in (7). 
Equation (7) means that the block inside of the big quark loop can be considered, with a good accuracy, as a block of 
the real process. This is a well-known feature of the high-energy jets; below we use it for estimating meson production 
amplitudes. 

We are not especially interested in the spectral integral (6) determined by the big quark loop, since it is quite common 
for the TT and p production. Our target is the spectral integral which corresponds to the amplitude T(j); qi, q2) that 
is the p or TT meson production block. This block is shown separately in Fig. 2b. 

B. Calculation of the loop diagram of Fig. 2b 

The amplitude of the loop diagram of Fig. 2b represented as a double dispersion integral is: 

M = T{p;qi,q2)2poi27rf6^''Hp + qi-q2-p'), (8) 



m/ \ f dsds f meson (•5) G^mesonf^ ) ri 

T[p;qi,q2) = / ^ / dcp 5 ; ^ Ameson • 

J TT J S — p, S — p-^ 

Here s and s' are the invariant masses squared in the intermediate gg-statcs; Gmcson is the vertex for the meson —^ qq 
transition. The ratio G'nieson(s)/(s — p"^) determines the wave function of the produced meson up to the spin factor, 
and S'mcson is the spin factor of the loop diagram 2b. 

Let us explain the structure of the spectral integral (||) in more detail. 

The cut quark loop should be integrated over the phase space of the intermediate qq-states, $(P;fci,fc2) and 
$(P'; k'l, k'2), with the invariant phase space determined as 

^f^^'-'^' - -2 i4fk; IX (^')**™(^ ' ■ (») 

The integration is carried out taking into account the conservation law for momenta which flow through the ladder 
(wavy lines in Fig. 2b); the momenta are denoted as qi and 92- Thus, the phase space integration of the cut diagram 
is given by the expression: 

d$(P; fci, fc2)rf$(P'; k[,k'2){2TTf2k[^5^^\]^\ - ki - qi)(27r)32fc^o5(3)(]^^ _ ^^0) 

Here kf — k'^ — and = s, P'^ = s' . We use light cone variables for the wave functions of the produced mesons; 
the simplest way to introduce them is to use the infinite momentum frame. Then the four-momenta P = (Pqj P±) Pz) 
and ki = {kiQ, kj_L , k^^) are written within the limits p — > 00, kiz 00, fc-^ 00. After introducing the hght-cone 
variables (^i = ki^/p, S,^ — k[^/p, mf^ = ni^ + k h mf^ = + kf^) and putting qi^ small {q'i^/p ^ 0, the constraint 
of the multiperipheral kinematics) we re-write (ttO) as follows: 

1 . ^1^5(1 - a - i2)d'k,^d'k2^5^^\\^^^ + k2^)5 (s-^ ^ 

167r ^1^2 V ii 4: 



?2 



^5[s' + (qi^ - ci2^f - ^ - ^ j • (2^)'2Po5(')(P + qi - q2 - P') 

= d^(a,6;kij.,k2i;qi^,q2±) • (2^)32Po5^''(P + qi - q2 - P') ■ (11) 
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The factor dif) stands for the phase space integration in the spectral integral (H): d(f> = ^2; ki_L, k2_L; qi_L, q2±)- 

In the approximation given by @), when the jet block inside of the big quark loop is considered as a real process, 
one may fix qi = q2 = 0; then the inclusive cross section is proportional to T{p; 0, 0). Taking (5-functions into account 



in (11), the formula for T(p; 0, 0) acquires a rather simple form: 
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Tb;0,0) = j^J / d^k^ ( ^—['^ ) 5„ , (12) 



167r3 J ^{1-0 J V 


where s = m^/ (^(1 — ^)). 

The amplitude r(p;0,0) alone does not determine the inclusive cross section da/dx{Z'^ — s- meson + X) because the 
amplitudes A{Wi,0) and A(W^|,0) depend on ^ and fc^, sec Taking into account this dependence, one has at 
a; ~ 0: 

^ / \ 2 



(Z" -> meson + X) ^ ^ J J d^k^ (^^™- ) 5„ n(l^f , Wi) . (13) 







The spin factor 5moson is closely related to the normalization of the wave function of produced meson. 

C. Spin factors Sp and 

Below, we calculate the spin factors Sp and S^r at qi = q2 = 0. They are: 

— — Sp {i^ii{kx + TO)(fc^ + m)i"f5{—k!2 + TO)(-fc2 + m)^ 
Sp= -Sp (7,^(^1 +m)(fc;+m)7^(-fc^+m)(-fc2+m)) . (14) 

We have taken into account here that the quark-gluon ladder carries quantum numbers of the scalar state, = 0+; 
hence the quark-ladder vertex is unity. 
The p meson vertex reads: 

la — flaa'la' i 9aa' — 9aa' • i-LOj 

In the spin factor Sp the summation is performed over polarizations of the meson. For the spin factors we have: 

S^ = Srn^s , Sp = lQ rr?{s + 2m^) . (16) 

Let us now demonstrate that similar spin factors determine the normalization of wave functions of the p-meson and 
the pion. 

D. Wave functions of the pion and the p- meson 

In the framework of the light-cone technique it is reasonable to introduce the wave function of a particle and its 
normalization using the form factor of the particle. The procedure of definition of the wave function is discussed in 
detail in |^ , ^ . Schematically, for the qq state this procedure looks as follows. 

The form factor of a composite system (for definiteness, we consider the pion form factor) is determined by the 
triangle diagram Fig. 2c, where the photon interacts with the composite system. The form factor is represented as 
a double spectral integral over masses of the incoming and outgoing pion; corresponding cuttings are shown by the 
dashed lines I and II in Fig. 2c. 

The structure of the amplitude of the triangle diagram for the pion has the following form: 

4*'-) = {p.+p',)F^{q^) , (17) 

where p and p' are momenta of the incoming and outgoing pions, the index v refers to the photon polarization and 
FT^{q^) is the pion form factor which, in terms of the double spectral representation, can be written as 
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J n-^ J s-fi^ 



4m ^ 



Here d(j)'^*'^\ki, k[.k2) is the phase space of the triangle diagram: 

(fci, fcl, ^2) = d$(P; ki,k2)d^P'; k^k'^) ■ {2nf2k'2^S^^^^ (k^ - ka) 



(18) 



(19) 



P and P' are total momenta of the qq system before and after the interaction with the photon: P = ki + k2 and 
P' — k[ + k'2, and P^ = s and P'^ = s' . The spin factor of the triangle diagram, TV, is determined by the following 
trace: 



(-)Sp [i-fsCK + ™)7^(fci + m)ij5{-k2 + m)J = {P^ + P'^^) T^{s, s' , q^) 
The index _L stands for vectors orthogonal to the photon momentum: 

At = one has ^^(0) = 1. Direct calculations of equation ( |l8|) in the limit ^ give: 



4m2 



(20) 



(21) 



(22) 



Here p{s) is the phase volume of the qq system: 



p{s) = - I d'^iP-MM) = 



and S^^'^ (s) is the trace of the quark loop diagram for the pion: 

S^:"f\s) = (-)Sp[z75(fci+™)z75(-fc2 + m) 

Using light-cone variables we come to the following form of (|2: 

1 

1 

1 = 



Am? 



2s 



d^k^L 



s — fl'' 



2s 



(23) 



(24) 



(25) 



where s = {n? + /s^)/ (^(1 — ^)). This equation enables us to introduce the pion wave function in the form 



■07r(C,kj 



s — p? 



(26) 



which is normalized by the standard requirement. 

Likewise, we introduce the p-meson wave function: it is defined by the form factor which is the spin matrix Faa' (ff')- 
In problems that do not deal with polarization properties of the vector particle, it is convenient to work with the 
trace of the form factor matrix, Faa{<f'), which is normalized by 



(27) 



a=l,2,3 



The trace '^aFaa{<f) is determined by the expression analogous to (^8|), with evident substitutions G^r Gp and 
Tp. As a result, we obtain the normalization for averaged form factor: 



CO 2 



(28) 



7 



where 



-is. 



la- Pa^j (fcl + m) i 7q - Pa^ 1 ("^2 + m) 



The p ^ qq vertex, 7^ — PaP/P^, selects three degrees of freedom of the p-nieson. We have 

Slr^'Hs) = lis + 2m') . 
In the infinite momentum frame, rt28) is re-written as: 



(29) 



(30) 



(31) 



where 



V'p(C,k_i 



(32) 



E. V/P ratio for the central region and colour degrees of freedom 

The normahzation conditions for pion and p-meson wave functions define unambigously the ratio of yields for prompt 
production: p/n — 3, if the wave functions of these mesons are similar. Indeed, the expression (|3|) represented in 
terms of wave functions ipp and V'tt give us (^ immediately. 

We have not taken into account explicitly the colour degrees of freedom in the derivation. This, however, can be 
easily done. For the meson qq vertex the colour operator is equal to I/^/l^, where / is a unity matrix in colour 
space. We have two colour amplitudes, singlet and octet, for the chain of the quark loop diagrams, Ai and Ag. The 
couplings of the amplitudes Ai and Ag to quarks {g{Ai) and g{As)) are proportional to / and A (Gell-Mann matrices). 
All colour operators are the same for both pion and p- meson production. Because of that, the colour factors are not 
important for the ratio p/n — they are identical and cancel in the production ratio. 

As was stated above, the main contribution into inclusive meson production comes from the ladder diagram As- 
This is due to the fact that the coupling constant for the amplitude with c = 8 is larger than for c = 1. In terms of 
1/Nc expansion g{Ai)/g{As) ^ I/VNc- 

III. INCLUSIVE PRODUCTION OF MESONS IN THE FRAGMENTATION REGION 



The inclusive production cross section of mesons in the fragmentation region is determined by the discontinuity 
of the diagram 2d. The spectral representation for this diagram is written as an integral over the masses of initial 
and final qq states in the transitions Z'^ —> qq and qq —>■ Z'^ and over qq masses in the transitions qq —* meson and 
meson qq. The amplitude of the diagram of Fig. 2d reads: 

00 00 
f dMHW gz{M')gz{M'') f dsds[ 

J 7r2 (M2-M|)(M'2-M|) J 7r2 ^ 



4m^ 4m ^ 

< dMkuk2,h)dMklk'2,K)A{W',{k2-k'^ 



meson 



y t^mcson V'-'/'-'mcson / 



(33) 



The vertices gz{M'') and gz{M''') are written for the Z° qq and qq Z° transitions. Spectral integrals over s 
and s' stand for qq — > meson and meson — s- qq (where meson means tt, p). The wave function ^/jmeson of the produced 
meson was introduced in an explicit form in Section 2 for the pion and the p-meson, and the factors d03(fci, k2, fcs) 
and d(f>3{k'i, k'2, fcg) define the integration over phase spaces in the left- and right-hand parts of the diagram of Fig. 
2d: 
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d(j)3{ki,k2,k3) 



2 (27r)32fcio (27r)32fc2o ' 

x(2.r^'^^(p-^.-^.) Ij^^i^^r^^'Hp-k.^k,). (34) 

Here ^ and ^ s. 

The block A (W^, (^2 — ^2)) defines the multiperipheral ladder (wavy line in Fig. 2d). This block depends on the 
momentum transfer squared (^2 — fcg)^ and the total energy squared W^: 

~ MKI-x) ; (35) 

X is the momentum fraction carried by the produced meson: x = 2p/Mz, where p is the longitudinal component of 

meson momentum, pmeson = (p + Ai5_/2p, 0,p). 

The spectra da{Z° meson + X)/dx fall rapidly with increasing x: this decrease is governed by A [W^, {k^ — k'2))- 
All the characteristics of (^^ listed above are the same for both pion and p-meson production, the wave functions 

( fr) ( fr) 

V'tt and ipp are also supposed to be the same. The difference may be contained in the spin factors 5^ and Sp 
which are 



Sif^'> = (-)Sp 



+ R-f5){k[ + m)i75(fc^ + m)(/c3 + m) 



X «75(fci + m)-f^{l + R-f5){-k2 + m){-k'2 + m) 



^y'-) = (-)Sp 



7^^(1 + i?75)(fcl + m)j'^^{k'3 + m){k3 + m) 



X 7^(fci + m)7j^(l + P75)(-fc2 + m)(-fc^ + m) 



(36) 



The factor 7^^(1+^75) is related to the vertex — > qq which is determined by the vector and axial- vector interactions 
(the ratio of coupling constants is ^ 2.63 for the u quark and ~ 1.43 for the d quark). For 3^^"^ the summation is 
performed over polarizations of the Z'^ boson (index v); 7^ is orthogonal to ki + ^2 and 'y'^ is orthogonal to fcj. + ^2- 
For the p-meson spin factor sj/^\ the summation is performed over polarizations of the boson (index ly) and 
p-meson (index a) with 7^ ± {ki + k^) and -f'^ ± (k[ + fcg) (see also ([l5|)). 

Normahzation factors Smeson{s) and Smcslnis') in ( p3| ) are related to the definition of meson wave functions, see 
( p6| ) and (^). The explicit expressions for spin factors (^ are presented in Appendix B. Using these expressions, 
we have an opportunity to perform numerical evaluation of ratios of prompt meson yields in the decay reactions with 
smaller energy release. In the case of the Z'^ decay, when Mz » m and Mz >> /i, a reasonably good approximation 
for ( p3| ) is to substitute the integrals over and M'^ by semi- residues in the poles = M| and A/'^ — M^. Then 
we have the following kinematics for real jets: 

ki^k[, k2 = k'^. (37) 
In this approximation we can put k^ — fcg and s — s' . Then we have for light mesons: 

oUr) oifr) 

-J^c.2Af|(l + P^), 'L^^GAfKl + P^), (38) 

O-jr Op 

which gives the ratio of the prompt yields p : tt = 3 : 1 in the fragmentation region x ^ 0.5 — 1 (more generally, 
V : P = 3 : 1 for hadronic decays Z*^ — > qq with q = u,d,s). 

The same ratio appears for the production of heavy quarks Z'^ = QQ, where Q — c,b. For example, in the case of 
h quark the spin factors are: 

2[Ml + 2ml + R\Ml^Aral)] , 

qifr) 



6[A4 + 2ml + R'^{Ml~Aml)] , (39) 



and thus the ratio Bp^^j^-i^^^ : Bprompt also equals 3 
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A. V/P ratio in the fragmentation region: comparison to data 

We have seen that V/P = i in the fragmentation region as well as in the central region. However, in the central 
region the comparison of quark combinatorics with experiment is hampered by the presence of a number of the decay 
products of highly excited resonances, while in the fragmentation region this contribution is suppressed by rapidly 
decreasing spectra. This means that the fragmentation region allows us to perform a model-independent verification 
of quark combinatorics. Here we carry out such a comparison based on the ALEPH data ||Tl| ]. 

To operate with meson spectra at a; ~ 0.2 — 0.8 we have fitted the spectra {l/atot)d<7/dx to the sum of exponents 
^CiC"'''^; the calculation results are presented in Fig. 3 for 7r^,7r'',p'' and {p,p)- The ratio of fitting curves drawn 
with calculation errors (shaded area) is shown in Fig. 4a for p/n. We see that for 0.6 < a; < 0.8 the data are in 
reasonable agreement with the prediction p/n — 3. 

The figure 4b, c demonstrates the ratios K*'^ / and K*^/K^: the data do not contradict the prediction, though 
the errors are too large to conclude anything more definite. 

IV. SUPPRESSION OF THE STRANGE AND HEAVY QUARK PRODUCTIONS 

In hadronic multiparticle production processes (in jet processes of the type of Z° —^ hadrons or in hadron-hadron 
collisions at high energies) the production of strange quarks is suppressed. Strong suppression is observed for the 
production of heavy quarks Q = c,b. One can guess that this suppression, being of the same nature for different 
reactions, is related to the mechanism of the production of new quarks at large separations of colour objects. This 
mechanism is seen in its pure form in the two-particle decays (the corresponding diagram is shown in Fig. 5a). The 
block of the production of a new qq pair in the two-particle decay is the same as that of meson production in jet 
processes (Fig. 5b), so it is reasonable to suppose that the suppression mechanism of the production of new quarks 
is similar for these processes. 

A. Two— meson decay of the qq state and soft hadronization 

The decay of the qq state takes place as follows: the quarks of the excited state leave the region where they were 
kept by the confinement barrier, and at a sufficiently large separation a new quark-antiquark pair will be produced 
inevitably: together with the incident quarks, these new quarks then form mesons (i.e. free particles). Schematically, 
this process (which is the leading one in terms of the expansion) is shown on the diagram of Fig. 5a: two quarks 
fly away (with the momenta pi and P2), and at large quark separations the gluonic field produces a new qq pair 
(the quarks with momenta k2, fca); then the primary quark (now with momentum ki) joins the newly-born one (fc2) 
creating a meson. Likewise, another newly-born quark (k^) joins the other primary quark (now with momentum k^) 
producing the second meson. 

The block with the quark-antiquark production, that is the transition 

qiPl) + q{P2) qih) + q{k2) + qik-i) + q{ki) (40) 

is the key process that determines the decay physics; it is shown separately in Fig. 5c. The process ( ^o|) is responsible 
for the leaving the confinement trap by quarks. In modelling the hadronization transition amplitude (pO|), quark 
combinatorics uses the hypothesis of soft hadronization. The idea was formulated at the 70's, and even now the soft 
hadronization hypothesis looks rather reliable and productive. It suggests that in the ladder of produced quarks (Figs, 
la, 5c) the contribution comes from small momentum transfers (of the order of R/^onfinement)- ^^"^ framework of 
the space-time picture this means that new qq pairs are produced at large separations, when colour objects leave the 
confinement well. 

The soft hadronization hypothesis applied to the decay processes treats the ladder diagram of Fig. 5c for the decay 
amplitude of Fig. 5a in the same way as for jet production, Fig. 5b: process ( ^ ) is an elementary subprocess both for 
the high-energy ladder and the two-particle decay amplitude, and the momentum transfers which enter the amplitude 
(|^ appear to be small in the hadronic scale, ~ R confinement ■ 

B. Spectral representation of the decay diagram 

Below we consider in detail the decay amplitude of Fig. 5a, performing calculations, as is done before, in the 
framework of the spectral representation with the light-cone wave functions for qq states. 



10 



By using the notations 

P =Pl+P2, ki2 =/si+/c2, ^34 =fc3 + fc4, 

M2 =(fci+fc2)', S12 =(fcl+fc2)', S34 ^{h + kif (41) 

we have the following spectral representation for the amplitude of the diagram shown in Fig. 5a: 

00 00 
^(99 state ^ two mesons) = / — — *„(M"')d$(P;pi,p2) / — 

X t{pi,p2; k2, k4)d^ki2; ki, k2) ^^'(^34; fcs, A:4)^'l(si2)^'2(s34) • (42) 

Here the masses of newly-born quarks i ~ 2,5 are denoted as ms, thus opening the way to consider the decay with 
strange quark production. The transition amplitude (|o|) of Fig. 5c is denoted as t{pi, p2; ki, k2, k^, k^). The decay 
amplitude (^) is written in terms of meson wave functions: for the initial state it is ^'m(Af^), and for outgoing 
mesons they are ^'1(512) and ^'2(512)- We do not specify here the meson spin structure of the quark propagator: we 
consider it below in a more detailed consideration of the transition amplitude fiO). 

Thus, the decay amplitude A is a convolution of the transition amplitude (MO) with wave functions of initial and 
outgoing mesons: 

A{qq state two mesons) = '^'i„ (g) t (g) ^'i\I'2 . (43) 
Further, we use light-cone variables; with pi^ — > 00 and fc^^ — > cx) we have 

Tn ~\~ k ? I 

ki — [kiz H , , kiz ) i — 1,4: 

TTl^ _j_ 

ki = {kiz H J J kiz ) i = 2,3 

'^Kiz 

2 



TTt -I— T)' 

Pl^{p^z + ^^^,P^±,P^z), i ^ 1,2. (44) 

2piz 

Let us use the frame where the outgoing particles are moving along the z-axis. We put 

P={Pz + ^^.Q.Pz). 

kl2 = ikl2z + 7T7^, 0, ki2z), 

2ki2z 

fc34 = (A:34. + :r^,0,fc34.). (45) 

The phase spaces (i$(fci2; ki, k2), d^lk^^; k^, k^) contain the energy-momentum conservation 5-functions which give 

^2 ^ m^+pli_ , m^+pl^ 



Sl2 



X\ X2 

w? + kl^ ml + k^j 



yi y2 



S34 = 11^^ + ^i. (46) 

2/3 Vi 

Here Xi = Piz/{.Piz + P2z) , Vi = hz/{kiz + fez) for i = l,2 and ?/i = hz/{ks.z + kiz) for i = 3, 4. 



C. The region of dominance of the transition amplitude tqq^qqqq 

Let us turn to the principal point, namely: the evaluation of the region in momentum space selected by the transition 
amplitude of Fig. 5c within the soft hadronization assumption. 

The hypothesis of soft hadronization means that the ladder diagram of Fig. 5c has a peripheral structure: it requires 
the momentum transfers squared to the qq block to be small, of the order of ^/ Rlonfinement ■ 
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-{Pl - fcl)^ ^ (Pl± - kl±)^ - ^confinement ' 

-(P2 - fc4)' ^ (P2± - k4i)' ^ ^;o'„/„„t • (47) 

Likewise, the momentum transfer squared in the quark propagator (see Fig. 5c) is of the order: 

{pi -ki- k2f ~ -(pii - kii - k2i)2. (48) 



Let us now discuss the peripheral constraint (|47|) in more detail. As an example, let us consider (pi — ki)'^. We 
have 

ip^-k.r^ip..-k,. + !!^-!!^)^, 

-(pii-ki^)2-(pi,-A:i,)2 . (49) 
The gluon carries a small fraction of the longitudinal momentum of the primary quark, therefore 

piz » Piz - kiz = A. (50) 

Then 

(fci -pif~ —{kf^ - pI^) - (ki^ - pii)2 ~ -(kx - Pi±)' . (51) 

Plz 

The peripheral constraint means that j/i >> j/2- But in numerical calculations this docs not require very large 
invariant energies squared, si2 and S34, which determine the blocks of quark fusion. For example, at j/2/yi ~ 1/10 
one has S12 ~ lOm^ ~ 1 GeV^ for m ~ 350 MeV, which gives for the quark relative momentum the value of the order 
of 300 - 400 MeV. 



D. The decay amplitude A(gg state two mesons) 

Let us now write the formula for i(pi,P2; fei, ^2, ^3, ^4) hi the simplest approximation, taking into account the 
t-channel propagators only: 

t{pi,P2;ki,k2,k3,k4) = ^ 



X - 



m2 + (pij_ - kij 
5^ (tos - 7(Pi-L - - k2j.)) 



+ (pi± - kii - k2_L)2 m2 + (p2J_ - k4^)2 



(52) 



To avoid ultrared divergence, the effective mass of the soft gluon is introduced into the gluon propagator (see, for 
example, [0). 

The equation ( |5^ ) does not state that the transition amplitude t{pi,p2; ki, /c2, fca, k^) selects large distances. At this 
point the amplitude (^) may be improved by incorporating form factors into the gluon emission vertex: 

5^5((Pi±-ki±)2), 5 ^ 5((P2± - k4±)2) . (53) 

With equations ( |5^ ) and (|3|) the amplitude A reads: 

1 

_ , f dXldX2S(l — Xl — X2) , , r/ X 

A(qq state two mesons) — / -—5 «Pi±ap2J.o(Pi_L + Pi±) 

J IQtT^XiX2 



J/l>J/2 



X / -—5 dki_LC(k2±(5(ki_L + k2_L) 

167r^?/i?;2 



X / — ^ dk3_Ldk4^(5(k3_L +k4_L)*m(xi,a;2;pij^,P2_L) 

g((Pi± -ki±)^) g^(m, -7(pi^ -kij. -k2±)) g((p2± - k4_L)^) 
ml + (pi_L - ki^)2 ml + (pij^ - kij_ - k2_L)2 + {^2± - ^i±Y 

x*i(yi,?/2;kij_,k2j_)^'2(j/3,?/4;k3j_,k4i) . (54) 
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This is the expression for the decay amplitude which makes it possible to discuss the rules of quark combinatorics. 

In a rough approximation that still gives a qualitatively correct answer, we neglect the momenta in the propagator 
of newly-born quarks: 

{nis - 7(pi_L - ki^ - k2_L)) ^ 
"ml + (pi_L - ki^ - k2j^)2 rris 

We have 

(y. 

A(qq state two mesons) = — • (g) t (g) ^'i^'2) , (56) 

rris 

This equation tells us that the probability to produce non-strange and strange quarks, uu : : ss = 1 : 1 : A, is 
determined by the ratio of masses squared of non-strange (w, d) to strange (s) quark. Introducing the constituent 
quark masses in the soft region, m„ ~ = m — 350 MeV and ~ 500 MeV, we get: 

9 

m 

A~ — -0.5. (57) 
m~ 

The equations (|5^) and ( p6| ) justify the statements of quark combinatorics applied to the decay processes [[7|-p|,p2|. Of 
course, we here suppose the identity of meson wave functions belonging to the same multiplet. 

The equation ( |57| ) gives us a rough evaluation of A, for in (^5|) we neglected momentum transfers squared that are 
compatible with light quark masses. In more sophisticated evaluations of A, one may take into account the momentum 
dependence of the quark propagator: 

(58) 



+ (Pi_L - ki^ - k2_L)2 ml+ < fc^ > ' 
where < /c^ > is a typical momentum squared inherent to the considered decay process. Therefore, 

A=4±^. (59) 

For standard decays of light resonances < >~ 0.1 — 0.3 (GeV/c)^, and this results in the increase of A compared 
to the estimation (p7|). Indeed, in the analysis A ^ 0.7 was found. Actually, equation ( pof ) demonstrates that A 
can vary depending on different types of reactions. 



E. Suppression parameter A in the hadronic decay 

Now let us turn to the processes of the qq pair production in jet processes Z*^ ^ qq —f hadrons, which is shown in 
Fig. 5b. All the above considerations, which have been used for the decay of a resonance into two mesons, may be 
applied to this process. As a result, we obtain the following formula which is a counterpart of Ipt 



A(Z° ^ two mesons + X) = — • (qqfromjet ladder g)^® \E'iVl'2) . (60) 



This formula differs from (|5^) by the initial state only, which is the wave function of qg-pair for jet ladder but not for 
the state defined by the wave function 5'i„. This means that the ratio of probabilities of producing a strange and a 
non-strange quark are given by the factor to^/to^. So, like in the decay process, one has: 
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m 



A~ — -0.5. (61) 



m 



Experimental data on the ratio of yields K^/ir^ do not contradict this evaluation. The ratio K^/ir^ as a function 
of X is shown in Fig. 4e. It is seen that at a; = 0.2 /n^ ~ 0.35. With the increase of x the ratio /n^ grows 
and reaches the value ~ 0.8 at x = 0.7. Such an increase is rather legible: the matter is that K mesons are produced 
both due to the formation of a new ss pair in the ladder (with the probability A) and fragmentation production of 
ss pair in the transition — > ss. Relative probabilities of prompt production of uu, dd, ss obey the ratio 

uu : dd : ss ~ 0.26 : 0.37 : 0.37. Because of that the production of K meson at large x is proportional to 
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A'+ - (0.37 • 1 + 0.26 • A) . (62) 

Let us comment this formula in a more detail. The meson is produced in the two processes: (i) fragmentation 
production of s quark (relative probability 0.37) and subsequent pick-up of the u quark from qq sea in a jet (relative 
probability 1), and (ii) fragmentation production of u quark (relative probability 0.26) and pick-up of s quark from 
the sea (relative probability A). 
The same quantity for pions is 

TT 

So, the ratio K^/n^ at large x is equal to 



(0.37 • 1 -I- 0.26 • 1) . (63) 



^^ 0.37 + 0.26A^ (64) 
7r+ 0.37-1-0.26 ^ ' 

at A = 0.5. This value agrees with the experimental data as it is demonstrated by Fig. 4e, where A(a;) determined 
as the i^^/7r^(x)-ratio is shown. 

The small value of /tt~^ at x = is a direct consequence of large probability to produce highly excited resonances: 
in the resonance decay more pions than kaons are produced. For the problem of breeding of strange and non-strange 
states in the decay, it is rather interesting to see the ratio K* / p — experimental data for K^* / are shown in Fig. 
4e too (shaded area). Remarkably, the ratio K^* / p^ has no tendency to decrease with decreasing x: this means that 
the rate of breeding of K^* and p^ in decays is approximately the same. Unfortunately, experimental errors are too 
large to have more definite conclusions about the behaviour of \{x). 

F. Production of heavy quarks 

The supression parameter for the production of a strange quark cannot be reliably determined. This is because 
the masses of strange and non-strange quarks are small compared to the mean transverse momenta of quarks in 
the production process, see (^9|). We can draw a more definite conclusion about the suppression parameter Aq for 
the production of heavy quarks Q — c,b. This parameter is defined by the same formula ( |5^ ) for multiperipheral 
production, so we have: 

Q 

Here we take into account that the gluon-quark coupling constant decreases with the growth of the quark mass, 
Aq ~ Q!^(mQ). The QCD scale constant, A, is of the order of 200 MeV. 

To estimate Ac and Ab, let us use m = 0.35 GeV, nic = Mj/^,/2 = 1.55 GeV, mj, = M-t/2 = 4.73 GeV and A = 0.2 
GeV. Then 

Ac~2.8•10-^ Afc ~ 1.1 • 10"". (66) 

The value of Ac should reveal itself in the inclusive production of J/vt and x mesons, while Af, is to be seen in 
reactions with T's: (^ J/^ -I- ^ x) + -'^ and ^ T -|- X. These reactions are determined by the processes 

^ cc^ c + {cc -I- gq— sea) + c and Z° ^ 66 ^ + (fefe -I- gg— sea) -I- b: for the production of a cc or hh meson a new 
pair of heavy quarks should be produced, since the quarks formed at the first stage of the decay, — > cc or Z^ — > 6&, 
have a rather big gap in the rapidity scale. So, within the definition 



Ac — r ycc c+ {cc + qq—sea) + cj /r(cc) 

Ab ~ F ^66 ^ 6 {bb + gg-sea) + b^ /r{bb), (67) 
we estimate r(cc — *■ c -I- (cc -I- qq— sea.) + c) and T{bb b + (bb + gg— sea) -I- b) by the available data from ||l^: 
r(^cc^c+ {cc + qq-sea) c^ ~ F (^J/^{1S)X) + T{J/-^{2S)X) + r{x{lP)X^ , 
F -I- (bb + qq-sea) -I- & ) - F ( T(1S')A: + T{2S)X + T{3S)x] . (68) 
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Experiment gives [lO[ |: 



T{cc)/T{hadrons) = 0.177 ±0.008, 
T{bb)/T{hadrons) = 0.217 ±0.001, 



r i^J/^'{lS)X + J/^{2S)X + x(lP)X j /T{hadrojis) = (1.17 ± 0.13) • 10"^ 

r (t{IS)X + T{2S)X + T{3S)X) /T{hadrons) = (1.4 ± 0.9) • 10"^ . (69) 



Thus, we have 



T\^J/'^{IS)X + J/*(2S')X + x(lP)Xj /r(cc) = (2.07 ± 0.23) • 10"^ 

rf T(1S')X + T{2S)X + T{3S)x) /r{bb) = (0.31 ± 0.19) • 10"^ , (70) 



in a reasonable agreement with (pq). 

V. THE PROBLEM OF SATURATION OF THE PRODUCED 
QQ AND QQQ STATES BY MESONS AND BARYONS; 
BARYON-MESON RATIO AND THE WATSON-MIGDAL FACTOR 

In quark combinatorics formulated in the baryon quark number reveals itself as the probability to produce a 
baryon containing this quark. For the hadronization of the quark qi in {q, q)-sea the production rule reads: 

q^ + {q, 9)sea ^ ^B, + ^M, + hi + (Af, B, B)sea, (71) 

where Bi and Mi are baryons and mesons containing the quark qi and M, B and B are mesons, baryons and 
antibaryons of the sea. 



Generally one can write 



M = Y,^lLML , B = Y.(3lMl , (72) 



L L 



where indices L = 0, 1, 2, ... define the multiplet, while /i^, //^'•' and /J^*'' are production probabilities of mesons 
and baryons of the given multiplet in the process of quark hadronization. These probabilities are determined by 
characteristic relative momenta of the fused quarks. 

A. The baryon-meson ratio 

Our present understanding of the multiperipheral ladder is not sufficient to re-analyse (|7^) on the level carried out 
in Sections 2 and 3 for V/P. Nevertheless, the data for decays Z° ^ p + X and Z° ^ 7r+ + X definitely confirm the 
equation (|7l|). Quark combinatorics [|| predict for p/7r+ at large x: 

p/7r+ ~ 0.20 (73) 

In Fig. 4d one can see the p/n'^ ratio given by the fit to the data (shaded area) and the prediction of quark 
combinatorics ([T^): the agreement at x > 0.2 is quite good. 

Let us comment the result of our calculation p/7r+ ~ 0.20 for leading particles in jets. In the jet created by a 
quark the leading hadrons are produced in proportions as it is given by eq. ([7l|): Bi : 2Mi : M. We consider only 
the production of hadrons belonging to the lowest (baryon and meson) multiplets, and, hence, keep in ([z^ only the 
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terms with L = (hadrons from the quark S'-wave multiplets). In our estimations we assume /3o — Mo 7 find therefore 
we substitute Bi B,{0), M,{0) and M M(0). The precise content of 5,(0), Mi{0) and M(0) depends 

on the proportions in which the sea quarks are produced. We assume flavour symmetry breaking for sea quarks, 
uu:dd:ss = 1:1: A, with < A < 1. For the sake of simphcity, we put first A = (actually the ratio p/ir depends 
weakly on A). Then for the u-quark initiated jet we have: 

B«(0) 

M„(0) 

M(0) 

The hadron content of the d-quark initiated jet is determined by isotopic conjugation p n, n p, tt"*" —f ■k~ , 
and the content of antiquark jets is governed by charge conjugation; in jets of strange quarks only sea mesons (M) 
contribute to the p/i:^ ratio. 

Taking into account the ratio Bi : 2Mi : = 1 : 2 : 1 and the probabilities for the production of quarks of different 
flavours g^, given by (1), we obtain p/7r+ ~ 0.21 for A = 0. We can easily get the p/7r+ ratio for an arbitrary A: the 
decomposition of the ensembles 5^(0), Mi(0), M(0) with respect to hadron states has been performed in Ref. (see 
Appendix D, Tables D.l and D.2). But, as was stressed above, this ratio is a weakly dependent function of A: at 
A = 1 we have p/n'^ ~ 0.20. 



2 1 .A 

—p + —n + (A — resonances ), 
15 15 

» ^7r+ + -;^7r° + T7:{''1 + 1') + { vector mesons ), 
8 Id Id 

— 7r+ + — 7r° + — TT" + T^iv + "n') + { vector mesons ). (74) 
Id Id Id Id 



B. Production of highly excited resonances and the Watson-Migdal factor 

For quark combinatorics the saturation of qq and qqq states by real hadrons is of principal importance. The 

probability of saturation is defined by the coefficients {fj,L, ^J■^l^) and (Pl^Pl^) in (72). The main question is what 
contribution from high multiplets is important enough not to be negligible in the spectra. 

Consider in more detail the production of mesons in the central region: qq — > M. The central production of qq 
states is provided by the diagrams of Fig. 6 (loop diagram. Fig. 6a, and interactions of the produced quarks, Fig. 
6b). The diagrams of the type of Fig. 6b for final state interactions lead to the relativistic Watson-Migdal factor. 
To estimate how many highly excited states are produced, we have to find out which states are determined by the qq 
system in the multiperipheral ladder. 

The constructive element of the ladder is a process shown in Fig. 5b. In this process, as was stressed in Section 
4, new qq pairs are created at relatively large separations (in hadronic scale), at r ~ 1 fm: these separations are just 
those in the qq —^ M transitions. The orbital momenta of the qq system for this transition can be written as i ~ kr. 
For relative quark momenta k < 0.6 GeV/c, we have: 

L < 3. (75) 

Relying on the behaviour of the Regge trajectory, one can understand, to what meson masses fi this relation corre- 
sponds. The trajectories for qq states are linear up to ~ 2.5 GeV Q: 

a{fi'^)~a{0)+a'{0)fi\ (76) 

the slope a'{0) is approximately equal a'(0) ~ 0.8 GeV^^ and the intercept belongs to the interval 0.25 < a(0) < 0.5. 
Hence, for large /x, the estimation gives ~ a{^'^)/a' (0); with a(M^) ~ 3, we have ~ 4 GeV^. So we conclude 
that in the multiperipheral ladder it would be natural to expect the production of qq mesons with masses 

< 2 GeV. (77) 



It is instructive to write down the qq states which belong to this interval and saturate the expansions ([7 21). The qq 
states are characterized by the orbital momentum L, the quark spin 5 = 0, 1 and the total momentum J; one more 
characteristics of the qq state is the radial quantum number n. In the region /-t < 2 GeV the following qq multiplets 
n'^^+^Lj with n = 1 are located [|lO|,||: 

(78) 



1^5*0 , 


l'5i ; 








l^Pj 


(J = 


0,1,2) 






(J = 


= 1,2,3) 






(J = 


2,3,4) 
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As ([75|) tells us, all these states contribute significantly to the meson production. 

The contribution of states with n > 1 is determined by the structure of the interaction of quarks in the final state, 
see Fig. 5b. For each partial wave the sum of diagrams with final state interaction is: 

l-i^,) <^\-B*{s'q,) ' (^^^ 
where 1/(1 — B{sqq)) is the relativistic Watson-Migdal factor: 

oo 

B{sqq)^ I ^ Re B{sqq) + iN{sqq)p{sqq) . (80) 

J n S - Sqq 
Am? 

Here the function N{s) characterizes the interaction of quarks for given partial wave. 

The cross section of qq pair production with invariant mass squared Sqq is determined by (^9|) at s'qq = Sqq: 

da(q-q{Sq,) niSq,, Sq,) ( ^ _ ^ (,^^_) ) ^ ^ _ 5. (,^^_) ) • (^D 

It follows from quark-hadron duality that da qq{sqq) + X ) describes in the average the spectrum of produced 



resonance: 



dcr yqqisqq) + Xj ~da\^'^2(Resonancesnearsqq)+Xj. (82) 

The averaging of the resonance production cross section is performed over a certain vicinity of Sqq. In this way we see 
that the structure of the Watson-Migdal factor determines the rate of decrease of da I qq{sqq) + X j with the growth 



of the invariant mass as well as the quantitative contribution of the resonances with n > 1. 
There are the following qq multiplets of the radial excitations in the region /i < 2 GeV ^] : 

n^So , n^Si (n = 2,3); (83) 
n^Pi , n'^Po.u2 {n = 2) . 

The probability to produce these states is determined by the decrease of da I qq{sqq)+X I , and by the Watson-Migdal 



factor in particular. 

Summarizing, we definitely expect a considerable production of resonances belonging to qq multiplets with n = 1, 
see (78). With less defiteness we can judge upon the production of resonances with n > 1 presented in (83): if the 
Watson-Migdal factor suppresses strongly the contribution of large Sqq, the production of multiplets is suppressed as 
well; but another situation is also plausible. 

Only a small part of resonances belonging to the multiplets (78) and located at 1500 — 2000 MeV is included into 
the tables of Q ; a considerable number of resonances were discovered and identified only recently, see [^3|-^ and 
references therein. Moreover, numerous decay channels related to multiparticle states have not been identified yet. 
This makes the restoration procedure of the prompt probabilities on the basis of exclusion the known resonances 
|^-|9) doubtful. 

There is one more effect which makes the realization of the program of ||37|-|39| questionable presently: the effect 
of the accumulation of widths of the overlapping resonances by one of them ]40| ]. An identification of broad states 
with widths of the order of r/2 ^ 400 — 600 MeV in the mass region 1500 ~ 2000 MeV looks rather ambiguous at 
the present level of the experimental data, although the existence of such a type of states seems to be probable if we 
expect the existence of exotic mesons (glueballs and hybrids) in this mass region (see Goi for details) . 



VI. COHERENT PRODUCTION OF THE cu AND p° MESONS AT LARGE X 

Indications on coherent phenomena in the hadronic Z*^ decays follow from the ratio w/p" at x ^ 1, see Fig. 4f. In 
multiparticle processes, where the yields of p and uj occur without interference of flavour components, the ratio 
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ujIp^ = 1 



(84) 



is valid. 

Because of the approximate equality of masses of p and w, the equation (|8^ ) is not violated by the decay of the 
highly excited resonances. The figure 4f demonstrates the validity of (Q) at a; ^ 0.1 — 0.2. However, at a; ~ 0.5 — 1.0 
the ratio is definitely less than unity, thus unambigously pointing to the interference of flavour components of 

the amplitudes uu and Z° dd. 

Consider the production of uj and p in the quark jet. If the hadronization of quarks in jets Z'^ u -\- X and 
— > d -\- X happens in a non-coherent way, the leading quark in the w-jet picks up an antiquark u from the 
sea, giving rise to the transition uu — > [uu + dd)/2 + [uu — dd)/2 — + p^ /^/2. Likewise, in the d-jet dd 

[uu + dd) /2 - [uu - dd) /2 = w/V2-p°/\/2. If there is no interference of u- and d- jets, the equation (^4|) is valid for 
each of them. But if the production amplitudes of uo and p do interfere, the eqaution ( p^ ) is violated. Let us consider 
the ratio lo/ p^ for such a case. 

The amplitude of fragmentation production of u and d quarks is determined by the vertex of the transition Z'^ —> qicji 
which has the following structure ['ipij^[gv — 5)i75)'0i] ^'^^ ^he u quark = 1/2, 5^ ' = 1/2 — 4/3sin^ 0^ and for 
the d quark g'-f = -1/2, g\f^ = -1/2 + 2/3 sin 9w Therefore, for the axial- vector interaction we have the following 
transition: 

A 0.5 [u + [qq - sea + u)] - 0.5 [d + [qq - sea + d)] (85) 

0.5 [[uu + Xuu] - 0.5 [[dd + X^^] . 
If Xuu and X^^ stand for the same state, Xuu — = X, then 

A^^[uu-dd)+X = V2oj + X . (86) 

The amplitude which is due to the vector interaction is 

V 0.19 [u + [qq - sea + u)] - 0.35 [d + [qq - sea + J)] (87) 



^ 0. 19 [[uu + Xuu]-0 -35 [[dd + X^^] . 

With Xuu = X^g we have: 

V V2[-0A6uj + 0.5Ap°) + X . (88) 
Therefore, for the axial-vector and vector interactions the ratios are: 

^) =0, (^) =8.8.10-. (89) 

" / axial V " / vector 

Hence, the coherent processes strongly suppress the production of to meson compared to p*^. In the experiment 
oj/p'^ ~ 0.4 at X ~ 0.5 — 1, that enables us to evaluate the contribution of coherent processes to be of the order of 
50%. 



VII. CONCLUSION 



The rules of quark combinatorics are the consequence of the quark structure of hadrons, which directly reveals 
itself in the ratio Pprompt/i^ prompt ~ 3. This equality is valid for the multiparticle production processes such as the 
hadronic Z^ decays and high-energy hadron collisions. 

The spectra observed in the multiparticle production processes are formed mainly due to the decay of highly excited 
resonances. In the meson sector states with masses up to 1500 — 2000 MeV contribute significantly; we can suppose 
that in the baryon sector it would be the contribution of states up to 2000 — 3000 MeV. The reconstruction of the 
spectra of highly excited resonances does not seem possible on the present level of knowledge. Therefore, one should 
investigate the prompt particle yields using jets with large x where the contribution of the resonance decay products 
is considerably suppressed. 
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The ratio P^/tt" = 3 observed in the hadronic decay at x '--^ 0.5 — 0.8 agrees with the predictions of quark 
combinatorics reasonably well. 

The value p/tt"*" is also in agreement with the quark combinatorial predictions at large x, thus revealing the small 
(about 25%) probability to meet the diquark loop in the chain of quark loops. 

The peripheral production of new qq pairs in the multiparticle production processes allows us to introduce the 
suppression parameter for the production of strange and heavy flavour quarks. The suppression parameters are 
determined by the ratio of quark masses, their values agree with the experimental data. 

In the central region, at a; < 0.2, the experiment gives us w/p" ~ 1 that also agrees with the quark combinatorial 
rules. However, at a; > 0.5 the experimental data definitely tell us that w/p" < 1 which is an indication to the 
interference of flavour amplitudes with the fragmentation production of u and d quarks: the signs of flavour amplitudes 
favour to / p*^ < 1 . 

As a conclusion, we can state that the predictions of quark combinatorics for hadronic decays of the Z° boson agree 
reasonably well with the experimental data. 
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APPENDIX A: V/P RATIO IN THE HIGH-ENERGY HADRON-HADRON COLLISIONS 

We demonstrate here that the ratio Vprompt / Pprompt = 3 is not specific for the decay processes — it is inherent in 
the hadronic multiple production processes as well. 

We consider the meson central production cross section, which is governed by the two-pomeron diagram as is shown 
in Fig. 7a. The probability of the meson production is determined by the structure of a pomeron ladder; here we 
use a QCD-motivated pomeron based on the pQCD pomeron of jl6j. In ||l^ the pomeron was constructed inserting 
the running coupling constant and soft interaction boundary condition. The pomeron found under these constraints, 
Lipatov's pomeron |16[| , has an infinite set of poles in the j-plane; it is a suitable object for using it as a guide in the 
analysis of the soft interaction region. Such an extension of Lipatov's pomeron into the soft region has been already 
used in |^ ; below we follow the results of . 

The gluon ladder which responds to Lipatov's pomeron is shown in Fig. 7b. For the description of soft interactions 
the quark loops can be incorporated into the gluon ladder since, according to the rules of 1/A^ expansion [|l^, the 
quark loops do not cause additional smallness. The cutting of quark loops (see Fig. 7c) provides the diagram which 
corresponds to the inclusive production cross section of the qq pair i.e. the production of meson states. In this sense 
the diagram of Fig. 7c is re-drawn as Fig. 7d: the quark-pomeron vertex is defined by the interaction of two reggeized 



gluons with a quark. The spin structure of such a vertex can be easily calculated (see 41 1, Appendix C) using the 
analysis of the leading-s terms in reggeon kinematics (large s and small t) for the vertex vector particle + fermion 

ii- 

In the leading terms of the 1 /N expansion there exist diagrams of the type of Fig. 7e in which the gluons of the 
upper (or lower) pomeron interacts with two quarks simultaneously (generally, these diagrams can be presented as 
Fig. 7f). As was shown in |^^, the sum of all these contributions (Figs. 7d, 7e etc.) results in the colour screening 
effect: the quark loop amplitude is nearly zero for the mentioned quark configurations at rgg ;^ 0.2 fm (the existence 
of the colour screening effects in hadron production processes was emphasized long ago [^^Q). 

The colour screening effect for the quark loop diagram in the pomeron ladder pl[ allows us to restrict the calculations 
of Vprompt/ Pprompt to the proccss of Fig. 7a, because this diagram is responsible for the main effect at large qq 
separations, rgg > 0.2 fm. 

The loop diagram of Fig. 7a is defined by the formula which is quite similar to (8) — only the spin factors are 

(P) (P) 

substituted as follows: ^ Sir and Sp ^ Sp , that is due to the accounting for vertices Pqg. The spin factors 



SP and read: 



Sl,^^ = -Sp (^ijsiki + m)h+{ki + m)i^z{-k2 + m)(-n_)(-fc2 + rn)^ 

SP = -Sp(7^(fci+m)n+(fci+m)7^(-fc2+m)(-n_)(-^2+m)) . (90) 
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Here n+ is the quark-pomeron vertex (the couphng to the upper pomeron in Fig. 7a) and (— n_) is the antiquark 
ponieron vertex (the couphng to the lower pomeron); we have also taken into account that pomerons in the diagram 
of Fig. 7a for the inclusive cross section carry zero momenta. 

If the upper pomeron in the graph of Fig. 7a moves along the 2;-axis with momentum pA = {po, 0,Pz) — {Pz, 0,Pz), 
while the lower one moves in opposite direction, we have 

n_ = (1,0,-1), n+ = (1,0,1). (91) 

As in Section 2, the loop diagram has been calculated in terms of the spectral integration; because of that k\ = k"^ = 
m^. After the operators h+ and n_ have been commutated with fci and k2, we have: 

5^^' = -Sp (i75(A;i + m)i^^{-k2 + m)) 4(n+fci)(n_A;2) 

Sf^ = -Sp (7^(fci + m)7^(-fc2 + m)) 4(n+fci)(n_fc2) (92) 

Up to a common factor, Si^^ and 5^^^ coincide with the normalizing spin factors of the pion and /O-meson wave 
functions given by (24) and (29): 

4^)=4(n+fci)(n_fc2)4"^) 

S<f) = 4(n+fci)(n_fc2) ^Sf'f^ . (93) 

This equation demonstrates that the two-pomeron diagram of Fig. 7a which defines the inclusive production cross 
section of p and tt in the central region of hadron-hadron collisions gives us Pprompt/T^prompt = 3. 



APPENDIX B: SPIN FACTORS FOR THE FRAGMENTATION PRODUCTION 

Below the explicit expressions for Sp^^ / S^^^ and Sy^^ / S^^^ are given for the processes ^ bb —^ B + X and 
Z° ^bb^ B* + X. 
For the pseudoscalar particle we have: 



Ai + i?2 ^2 , (94) 



where 



Ai = ((mt - m)2 - s') (M^ + 2ml) + (("^fc - - s) {M'^ + 2ml) , 

,2 \ / ^2 



— 3m^ ((mft — m)^ — s') — 3m^ ((mb — m)^ — s) ; (95) 
and for the fragmentation production of the vector particle: 



hr, = T77^3 + — r^4 + — r(^3 + A^) 



C>."'/^ SS SS SS' 



- {M^{M'^ - ml)Az + M'%M^ - ml)Ai - 3M'^M^ml{A3 + A^)) . (96) 



Here 



A3 = (s + s'){ml - TO^)^ + ss'(3m^ -ml - lOmbm - 4s') 

+ s'^{3ml — 2TOfcTO — TO^) , 
^4 = (s + s')(m^ — TO^)^ + ss'(3m^ — — lOmhm — 4s) 

+ s^(3m^ - 2mbm -w?) . (97) 
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FIG. 1. Hadronic decay of Z''^ meson, a) Multipcripheral ladder which provides the transfer of colour from quark to antiqurk. 
b) Self-energy diagram; the cutting along hadronic states (dashed line) determines the hadronic cross section, c) Diagram for 
the inclusive meson production cross section in the central region, d(7{Z° — > meson + X)/dx at a; ~ 0. d) Diagram for the 
inclusive meson production cross section in the fragmentation region at a; ~ 0.2 — 1. 
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FIG. 2. a,b) Production blocks for meson in the central region; the dashed lines show the cuttings of diagrams in the spectral 
integral, c) Form factor diagram determining the meson wave function, d) The cuttings of the diagram for the inclusive meson 
production cross section in the fragmentation region. 
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FIG. 3. The spectra {l/ahadron) dazo^meson+x/dx [11] for 
functions SOj exp(—6ia:) [11]. 
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FIG. 4. Ratios p°/7r" (a), K*^ / (b), K*^/K^ (c), p/7r+ (d), K*'' / and ;^±/7r± (c), Lo/p° (f) obtained from the 
exponential fit to ALEPH data [11] (shaded areas). Solid lines stand for the predictions of quark combinatorics rules. In Figs. 
4d, 4e, 4f the ratios p/tt"*", JsT^/tt^ and a;/p° which are obtained using a histogram description of the spectra [11] are also 
shown. 
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FIG. 6. Diagrams for the central production of qq pair. The set of loop diagrams responsible for the Watson-Migdal factor. 
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FIG. 7. Diagrams for meson production in the central region in high energy hadron-hadron collisions. 
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